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Abstract 



> 

, We study the quantization of many-body systems in two dimensions in rotating coordinate frames 

using a gauge invariant formulation of the dynamics. We consider reference frames defined by linear and 
quadratic gauge conditions. In both cases we discuss their Gribov ambiguities and commutator algebra. 
We construct the momentum operators, inner-product and Hamiltonian in both types of gauges, for 
systems with and without translation invariance. The analogy with the quantization of QED in non- 
covariant gauges is emphasized. Our results are applied to quasi-rigid systems in the Eckart frame. 



Oh! 1 Introduction 

-i— > ■ 

The problem of quantizing a many-body mechanical system in a rotating reference frame is of interest both 
by itself and for its possible applications to specific problems in, e.g., molecular and nuclear physics. When 
the underlying dynamics are rotationally symmetric, which is the only case we consider, the coordinate 
transformation from a space-fixed reference frame to a rotating one with the same origin is a time dependent 
symmetry transformation. It is thus appropriate to formulate the theory in such a way that it is invariant 
under symmetry transformations whose parameters depend on time, or gauge transformations pp. If the 
dynamics are described in terms of a gauge-invariant action, since we know how to quantize a mechanical 
system in a space-fixed coordinate frame, we can perform a gauge transformation in order to obtain the 
quantum theory in a rotating reference frame. Gauge invariance guarantees that both theories are physically 
equivalent. 

Whereas formulating a quantum theory in a rotating frame is a mathematical problem, the choice of the 
particular frame in which to formulate the theory is dictated by the physics of the specific system under 
consideration. It is often the case that the relevant rotating frame is defined implicitly, by restrictions on 
the trajectories of the system in that frame. In the gauge-invariant approach to the quantization in rotating 
frames, such restrictions are incorporated into the theory as gauge conditions. The action is then given in 
terms of degrees of freedom that are not independent, but must satisfy certain functional relations. This 
situation is familiar from the theory of gauge fields In quantum electrodynamics (QED), for instance, 
the degrees of freedom are the components of the vector potential A(t, x), which may be required to satisfy 
such relations as V • A = (Coulomb gauge), or x ■ A = (multipolar gauge) 0, at all times t. 

In this paper we study the quantization of many-body systems in two dimensions in rotating coordinate 
frames, using a gauge-invariant formulation. We consider systems of N spinless particles in the plane, 
interacting through two-body central potentials. We focus on developing the formalism, which is a necessary 
previous step to considering applications to realistic models. By restricting ourselves to the two-dimensional 
rotation group we separate the treatment of the gauge-invariant formalism from the technical intrincacies of 
non-abelian groups such as the three-dimensional rotation group, which we will consider elsewhere. 

* E-mail: abouzas@mda.cinvcstav.mx 
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Our treatment closely follows the approach of |21 E] to Yang-Mills theories. Previous treatments of the 
quantization of two- and three-dimensional TV-body systems in rotating frames within a gauge invariant 
approach have been given in El El E| (and references therein) . In 6 the gauge symmetry is implemented 
within the Hamiltonian BRST formalism |5]. In [7] a point of view based on the shape-space theory of 
deformable bodies is adopted. Non-gauge-invariant treatments can be found in, e.g., |9| in the context of 
nuclear physics, and in in molecular physics. Our goals are to establish a formal framework with the 
most direct geometrical and physical interpretation, and to provide careful and systematic derivations of our 
results. Furthermore, the results we obtain are different from previous ones. 

In the following section we describe the class of models considered throughout the paper, and their gauge- 
invariant formulation. Their quantization in a space-fixed frame is given, and shown to be equivalent to the 
original, non-gauge-invariant system. In the remaining sections we obtain the quantized theory in rotating 
frames by means of a gauge transformation from the space-fixed frame. In section El we consider rotating 
frames defined by linear gauge conditions. Such gauges are the most important ones for practical applications. 
We discuss the Gribov ambiguities of those gauges, which are important for the construction of the inner 
product in Hilbert space. The algebra of commutators is discussed in detail, and an explicit realization of 
that algebra in terms of differential operators is given. We then use those operators to construct the 
Hamiltonian. The equivalence of the quantum theory in these and other gauges with the theory formulated 
in a non-rotating reference frame is kept manifest at every step. We emphasize that by describing the 
system from a rotating reference frame defined by imposing restrictions on the coordinates, we are in fact 
introducing orthogonal curvilinear coordinates in configuration space. 

In section 01 we go through the same steps, though more briefly, to obtain the theory in the instanta- 
neous principal axes frame as an example of a gauge condition depending quadratically on the coordinates. 
Although the treatment is straightforward, the results are technically much more complicated than in the 
linear case. This fact makes the practical usefulness of this gauge condition doubtful. The quantization in 
linear gauges from section El is extended in section El to translationally invariant systems in rotating frames 
with origin at the center of mass. These results are then applied to the case of quasi rigid systems in section 
El where we discuss the Eckart gauge and recover some of the classic results of J21- I n section we give our 
final remarks. We try throughout the paper to make manifest the parallel between our approach and the 
quantization of QED in non-covariant gauges. In appendix El we give a brief summary of those aspects of 
QED which are relevant to the analysis presented in the main body of the paper. In appendix iBl we derive 
some technical results needed in the sequel. 

2 iV-particle system 

We consider a system of N particles in two dimensions interacting through a two-body central potential, 
described by the Lagrangian 

N 

£f(r Q ). (1) 

If the one-body potential U = 0, C is invariant under the group of Euclidean motions of the plane. In this 
and the following sections we consider U ^ and focus on the abelian group of two-dimensional rotations, 
deferring the discussion of translation invariance until section El 

We adopt the passive point of view for coordinate transformations. C is invariant under time-independent 
rotations of the coordinate frame. In order to make C invariant under changes of arbitrarily rotating 
coordinate frames we apply the usual Yang-Mills construction [13 to QJ. We add a new degree of freedom £ 
to the system, and postulate the following transformation law under infinitesimal rotations of the coordinate 
frame, 

5r a = -59zAr a , 5£ = -56, (2) 

with 56 = 56{t) an arbitrary function of t and z a unit vector orthogonal to the plane. These are the 
infinitesimal gauge transformations of the system. We define the covariant derivative D t r a = r a — £z A r a , 
which transforms as a vector under gauge transformations, 5(D t r a ) = —59z A D t r a . Substituting time 



^ N N 

a=l a<3=l 
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derivatives in by covariant derivatives, we obtain a Lagrangian invariant under time-dependent rotations 
of the coordinate frame. Explicitly, we write, 

1 N 1 N . 2 N N 

£ = 2 m »(D t r a ) 2 - V + 4C = g S maf "« + V mQr « ~ ^ ' X) m ° ( rQ A f ") " V + ( 3 ) 

Q — 1 a — 1 a — 1 a—1 

where we denoted by V the potential energy for brevity. In we added an extra term i z ^ to L which fixes 
the value of the angular momentum through the equation of motion for £. That term plays a role analogous 
to the source term j (x)A°(t, x) in electrodynamics 1 , as discussed in appendix^] L is invariant under gauge 
transformations if the constant £ z = 0, and quasi-invariant otherwise, 6C — —£ z 59. 

C in © describes the same dynamics as , but from a coordinate frame rotating with angular velocity 
— £ with respect to the laboratory frame ^3). Notice, however, that £ is a dynamical variable describing the 
coupling of the particles to inertial forces. The equations of motion for r a are m a D t D t r a + V a V = or, 
more explicitly, 

mprp = 2mp£,z A rp + mp£,z A rp + mp£ 2 z A (rp A z) — VpV, (4) 

where the terms corresponding to the Coriolis, azimuthal and centrifugal forces are apparent |15| . A conse- 
quence of rotational symmetry is the conservation of the system's total angular momentum 

N 

L z = z-^2 m a(r a A D t r a ). 

The equation of motion for £ is then L z — l z = 0. 

Since the system is gauge invariant we can fix the gauge by imposing a condition of the form 25({r Q }, £) = 
0, which is equivalent to selecting a rotating frame in which the trajectory of the system ({r a (i)},£(t)) in 
configuration space is constrained to satisfy the relation (5({r Q (t)},£(i)) = 0. The function © can be 
chosen arbitrarily, as long as any trajectory ({r' a }, £') can be transformed into a new one ({r a }, £) satisfying 
(3 = 0. The new trajectory must be unique, in the sense that no other trajectory obtained from ({r' a },£') 
by a gauge transformation satisfies the gauge condition. Otherwise the gauge is said to be ambiguous 
Supplementary conditions must then be imposed to fix the ambiguity. 



2.1 The laboratory frame 

Given any trajectory of the system ({r a (t)},£(t)) by means of a finite gauge transformation r' a = U(0(t))r a , 
£'(£) = £(t) — with U an orthogonal 2x2 matrix and 6(t) = j to dr£(r), we can obtain a physically 
equivalent trajectory with £'(t) = 0. The gauge condition £ = corresponds to choosing a non-rotating 
reference frame, the laboratory frame. In this gauge the Lagrangian © reduces to 0). The equation of 
motion for £, 53a=i ra a r a A r a — l z =0, which cannot be obtained from p|. must be imposed on the system 
as a constraint [2] . In the Hamiltonian formulation in this gauge that constraint is first-class |lf>| , not leading 
to further secondary constraints. 

The quantization in the gauge £ = is canonical. In units such that h — 1 we have, 

N N 

n = J2^—Pl + v > [rca,Pfij]=iS a f,5i j , (</#)= / Yl^r^* ({r a })<tp ({r a }) , (5) 

a=l TOq J (3=1 

with p a — l/iV a . The first-class constraint is imposed on the state space |16| . L z \ip) = £ z \ip). The constant 
l z can only take integer values in the quantum theory. We see that both the classical and quantum theories 
for this model in the gauge £ = are completely analogous to electrodynamics in Weyl gauge |2 El ( see 
appendixEJ. The constraint fixing the value of L z , in particular, plays the same role as Gauss law in QED. 

1 We remark that external source terms break gauge invariance when the gauge symmetry is non-abelian. In that case, in 
order to preserve gauge invariance and, at the same time, to have a non-homogeneous Gauss law, we must incorporate the 
source into the theory as a dynamical degree of freedom. 
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3 Linear gauge conditions 



The simplest form of gauge condition involving the coordinates is a linear relation among them. As discussed 
in section|S|this kind of gauge condition is relevant in the context of perturbative or semiclassical expansions. 
The following notations will be used throughout this paper, 



JY 



e ({ r a}) = ^2 m f3 {Apxp + B p yp) , (6a) 
0=1 

N N 

&{{r a }) = ^m a {B a x a -A a y a ), SH 2 = ^ m a (A 2 a + B*) . (6b) 

a—l a—1 

The general linear gauge condition is of the form 6=0, with $H 2 ^ 0. We denote position vectors R a and 
their components X a , Y a in this gauge by capital letters, as opposed to vectors in the gauge £ = (the 
laboratory frame) denoted by r a . Thus &({R a }) = but, in general, 6({r a }) ^ 0. This gauge condition 
selects a reference frame rotating in such a way that the linear combination of coordinates 6 vanishes for 
all t. If we choose, for instance, all coefficients in l|6a[) vanishing except for B%, the coordinate frame must 
rotate together with particle 1 so that particle stays on the X axis for all t. The formalism in these linear 
gauges is entirely analogous to that of electrodynamics in Coulomb gauge (appendix 0, in which the fields 
are also constrained by a linear relation. 

The transformation from the gauge £ = to the gauge 6 = is given by, 

R a (t) = U(6(t))r a (t), £(t) = -0(t), with 6(t) = arctan + mr, (7) 

where we denoted s = &({r a }), q = £2({r Q }) for brevity. The indetermination in 9 up to addition of ir 
is a Gribov ambiguity (see also |21 E]), related to the two possible choices £2 = ±\/q 2 + s 2 . We fix 
the ambiguity by requiring Q > and — tt < 8 < ir. Due to the relation £ = —6 for all t in this gauge, 
with 6 from Q, we can use {R a ,0} as dynamical variables instead of {R a ,£,}- The former set of variables 
is preferable in the operator approach we espouse in this paper, whereas in the path integral formulation 
switching from one set to the other amounts to a mere change of integration variables. The inverse to the 
transformation (JTJ is then, 

r a (t) = U(-e(t))R a (t), (8) 
with 9 an independent variable. (If we choose £ as a dynamical variable instead of 9, then in (JHJ) 9(t) = 

-JL*rt(r).) 

Solving the equation of motion for £ we get, 

/ N \ 

A R a - 4 J • (9) 

Since the gauge has already been completely fixed, we can substitute JSJ) back into the Lagrangian Q). 
Notice that the vector product appearing in is not the total angular momentum of the system. In the 
classical theory we obtain the momenta Tl a conjugate to R a by differentiating C in @ with respect to R a 
under the constraint 6({R a }) — &({R a }) = 0, to obtain, 

IIx a = m a X a + m a £_ ^Y a + > n Y a = m a Y a - m Q £ ^X a - ~^r^j > ( 10 ) 

where £ is given by 10. These momenta are consistent with the gauge condition, since they satisfy 




^ lip = 6 ({n a /m a }) , (11) 
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the last equality following from the linearity of 6. Relation is analogous to the condition that the 
momentum conjugate to the potential in Coulomb gauge must be transverse, eq. (|A.12|I . 

From Q and (JSJ) we can obtain the relation between the velocities {r a } hi the gauge £ = and those 
in the gauge 6 = 0, {R a ,9}. Correspondingly, we can express the momenta {p a } in one gauge in terms 
of momenta {II a } and l z in the other. With those transformations we obtain from H in JSJ) the classical 
Hamiltonian in this gauge, 

N 1 5K 2 N 

w = E^r iT ' + ^(^- A ) 2+v ' with A = J2( x n ^- Y ? n ^)- ( 12 ) 

a=l Q /3=1 

The quantity A defined by this equation will be henceforth referred to as the "residual angular momentum." 

By the same token, expressing TI a and R a in terms of p a and r a and using the Poisson brackets JSJ we 
get the Poisson brackets in this gauge. Alternatively, they can be found as Dirac brackets relative to the 
second-class constraints &({R a }) = = 6({i7 Q /m Q }). The result is, written in the notation of quantum 
commutators, 



(13) 



All other commutators among components of R a and iT^ vanish. The correspondence between (|13|l and 
(|A.14(I is apparent. Using (|13|) we obtain the commutators for A, 

[X a ,A] = -iY a -i^, [Y a ,A]=iX a -i^, 

m * (14) 

[A, iTxJ = iII Yj + i^^Q ({n a /m a }) , [A, 77k,] = -*I7x, + i—^^O. {{n a /m a }) . 

Furthermore, from 113|) and (|14|l we get, 

[6, = [6, = = [6, A] = [6 ({i7 a /m Q }) , A], 
[0,2T X J =im a B a , [Q, n Y J = -im a A a , [£J, 6 ({JT Q /m Q })] = 0, [0,A] = -i6. 

We see from (|13pi — (|15|) that JT Q and A generate translations and rotations, respectively, of {R a } on the 
surface 6 = 0. We see also that the gauge condition 6 = and ((TTJl are operator equations, that can be 
evaluated within commutators. 

In the quantum theory a realization of the commutators (|13|l is obtained by defining II a as the projection 
of the gradient V Q on the hyperplane tangent to the surface 6 = (which in this case is the surface itself, 
since 6 is a linear function), 



_ 1 1 m a A a A / 
Ux - - 7 dX a 7^~ lj \ A "dX p + ~dYf) 

_ 1 1 m a B a A / 



(16) 



These operators, which are analogous to (|A.15(I in QED in Coulomb gauge, satisfy the equation 6({II a /m a }) 
0. From (|12[1 and (|16|) we obtain the expression for the residual angular momentum operator A, 



f}=l 



which is shown in appendix [B] to have integer eigenvalues. 
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Using the relations Q and © between {r a } and {R a ,0} and applying the chain rule we obtain, after 
appropriately rearranging the derivative operators, 

Pxa = -— = cose (n Xa + -^( L * - A ) j - sin # [ n Ya + -^( L * - A ) j » 

Pya = 1 A = sinfl (i7x Q + ^^(L, - A)) + cos0 (lI Ya + ^^{L z A)) , 

with L z = 1/id/de. The first-class constraint L z ip = £ z ip is trivial to solve in this gauge, ip({R a },0) = 
ip({R a }) exp(i£ z e)/\/2ir. An analog of fBfy in QED is the simpler relation l|A.18|l . 

The Hamiltonian operator in this gauge is obtained from Ji in the gauge £ = as given by Q, through 
the transformation rules © and l(T%|) , 

« = g ^ (g^filfr, + 5^,^,) + f^(4 - A)» + V. (19) 

The first term in 7i has a structure similar to that of the Laplacian in curvilinear coordinates, with as the 
Jacobian and 77 as derivative operators. (In the case N — 1 the similarity turns, in fact, into an identity, 
see below.) Notice, however, that we did not postulate (JT5J, rather, we derived it from the expression JSJ) in 
the lab frame. 

The inner product in Hilbert space can be found from the expression Q for (4>\ip) in the gauge £ = by 
the familiar Faddeev-Popov method ^B], the relevant resolution of the identity being in this case, 

1= / daS(e({U(a) rfl }))e(£l({U(a)r p }))£}({U(a)r^), (20) 

J — 7T 

with 0(0) a step function enforcing positivity of the Faddeev-Popov determinant 0. As mentioned above, 
the condition > guarantees that there is only one root a = e (with e from Q), and not e + tt, to the 
equation 6 = in H20|) . Inserting i|20|) in the expression (J5J for (4>\ip) we get, 

(M>> = / II d 2 R 8{&)® {&)&<!>* {{R a }W{{R a }), (21) 

^ /3=1 

where we dropped a factor of 27r, the measure of the group SO(2). The hermitianity of A and V with 
respect to the inner product (|21|) is immediate in view of the commutation relations. In order to check the 
hermitianity of the Hamiltonian i|19f) it is then enough to show that the first term in (|19f) is hermitian. That 
calculation is straightforward, though somewhat lengthy, so we omit the details. 

It is sometimes convenient to redefine the state space by absorbing the Jacobian in the wave functions 
and eliminating it from the integration measure in the inner product. The redefined wave functions are 
if) = Q 1 / 2 ^, leading to the Hamiltonian, 

H = Q^HQ-V* = g ^- (n*, + n$ fi ) + g(4 - A)» + V - J*, (22) 

the last term being the quantum-mechanical potential [21 E]- Equation J22J), with classical II and A, is 
the Hamiltonian found in this gauge in the path-integral approach. Since the transformation Q depends 
non-linearly on {r a }, the associated change of integration variables in the generating functional entails a 
change in its discretization |19j . which ultimately gives rise to the quantum potential term. 

3.1 The case N = 1 

The case N = 1 is instructive EJ ■ By means of a time-independent rotation we can always reduce the 
condition & — to Y = 0. We thus fix a reference frame rotating together with the particle, so that it is on 
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the X axis for all t, with X > 0. In order to specify the position of the particle, we give its coordinate X 
and the angle 9 of the X axis relative to the laboratory x axis. We are then describing the motion in terms 
of polar coordinates with X the radial coordinate. 

From HI 6[l we have IIx = —id/dX and Hy = = A. The Faddeev-Popov determinant in this case is 
= X, and the Hamiltonian l|19(l reduces to that of a particle in polar coordinates, with angular momentum 
£ z . Similarly, the inner product i121ll corresponds to polar coordinates. Due to the constraint L x ip = £ z ip the 
integration over the angle variable 9 is trivial, so only the radial wave function appears in 121JI . If the wave 
function in Cartesian coordinates is ip(x,y), the radial wave function is ip(X,Q). The quantum potential 
term in also reduces in this case to its well-known form [312] f° r polar coordinates, — l/(8mX 2 ). 



4 Quadratic gauge condition: the instantaneous principal axes 

The quantization of the system in a rotating frame defined by a quadratic gauge condition follows the 
same lines as the linear case studied in section [3J Both the treatment and the results are technically more 
involved, however, so instead of considering a general quadratic gauge condition we restrict ourselves to 
the particular case of the instantaneous principal axes frame. That reference frame plays a central role 
in the treatment of rigid body dynamics. In the case of many-body systems, their quantization in the 
instantaneous principal axes frame has been proposed as a method for separating the "collective" rotations 
from the "intrinsic" dynamics. We briefly discuss that issue in section In this section we compute the 
quantum Hamiltonian and inner product by means of a gauge transformation from the gauge £ = 0. 
We define the quantities, 

^ N N N 

a— 1 a— 1 a — 1 

R 2 is the trace of the inertia tensor of the system, whose traceless part is given by (Z-'s qJ- The instan- 
taneous principal axes frame is then defined by the condition S = 0. As above, we denote vectors referred 
to this frame with capital letters, so that S({R a }) = 0. The gauge transformation from the gauge £ = to 
the gauge 5 = has the form (7J, with the parameter, 

m =^(&) + h (24) 

Due to the fact that the inertia tensor is second rank the number of Gribov ambiguities doubles with respect 
to the linear case {7J, there being now four solutions in the range — tt < 9 < tt. We fix the ambiguity by 
requiring Q({R a }) > and < 9 < tt. 

The Lagrangian is now (JHJ with £ having the same form as in . Consistency with the gauge condition 
S({R a }) — requires the velocities and momenta to satisfy, 

N 



m a {X a Y a + Y a X a ) = 0, (25a) 

Q = l 

N 

J2 (x*n Ya + Y a n Xa ) = o. (25b) 



We obtain the conjugate momenta in terms of velocities by deriving C with respect to R a under the constraint 
(|2"5aJ) . 

n x . - ™.i + {m „y„ (i + Mlgi) , n,„ = mA - w„ (i - f|£2) . (26) 

The basic Poisson brackets in this gauge can be obtained as in section |3J Since the momenta II a generate 
translations on the curved hypersurface S({R a }) = 0, which do not commute, the Poisson brackets among 
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momenta do not vanish. Correspondingly, in the quantum theory the operators iT Q do not commute with 
each other. The non-vanishing commutators among coordinates and momenta are, 

=-i%^, [Y ,n Xy ]=-i^L, (27) 

- 1 



[n Xp ,n Y ^] = ^ (m 7 x 7 n Y? - mpYpiix^) , [n Yfil n Y ,] = (m 7 x y n Xl3 - m x ? n x ^ 



Like in the previous section, we obtain a realization of this commutator algebra in terms of differential 
operators by projecting the gradient operator on the hyperplane tangent to S = 0, 

19 1 m a Y a A / d d \ 

nx ° -lax:- l^T fj + ' 
i a i m a x tt / a 9 x 

^ = - t-^- ^ + 

Both the classical momenta (|26H and the quantum operators l|28[> satisfy the constraint (|25b|) . The gauge 
condition S({R a }) — and its counterpart (|25b(l are operator equations, that can be evaluated inside 
commutators as can be easily checked from ()27|) . The Hamiltonian operator is written in terms of momentum 
operators (|28[) as, 



with R = V-R 2 and 

a . ± (x.n Y . - W . (i - *) £ - (i + f ) ± rjjfc (30) 

a— 1 v 7 a— 1 v 7 a— 1 

the residual angular momentum in this gauge. From l|29|l we find the form of the quantum potential, 

R 2 7 -AN 



Q 8Q 2 8i? 2 ' 

In order to find the inner product in this gauge we start from the resolution of the identity, 



(31) 



1=/ daS(S({U(a)r p }))Q(Q({U(a)r }))2Q({U(a)r p }). (32) 
Jo 

Notice that we restricted the integration range to < a < n. Alternatively, we can integrate from — n to 7r, 
and set the l.h.s. of l|32|l equal to 2. The inner product is then obtained as, 

N 

n 

7 (»> 



(<t>\1>) = \ / I] d 2 R p S(S)e(Q)2Q ( j ) *({R a })^({R a }) 



13=1 



where we divided by 2-7T, and omitted the argument {R a } in 5, Q and R for simplicity. In the second line 
we wrote the gauge condition as S/R, taking into account that the factor 1/R has no finite zeros and that 
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its singularity at the origin is suppresed by the zeros of S and Q there. This last expression for ((f>\ip) is 
convenient to check the hermitianity of 7i. We again omit the details of that proof. 

As a consequence of the gauge condition S({R a }) = being quadratic, the momentum operators Tl a 
in l|28|) have coefficients which are ratios of quadratic polynomials in the coordinates, whereas the operators 
(|f 6 II for a linear gauge have constant coefficients. Accordingly, the basic commutators (|27|l are rational 
functions of coordinates and, furthermore, momentum operators do not commute with each other. The 
structure of the Hamiltonian operator in this gauge, eq. (|29|) . is also much more complicated than in linear 
gauges, eq. (I19II . As shown in appendix IE! the residual angular momentum operator A, l|3()[l. does not have 
integer eigenvalues. 



5 Center of mass motion 

In this section and the next one we set U — in the Lagrangian and take into account the translation 
invariance of in order to separate the center of mass degrees of freedom. Since the motion of the center of 
mass is dynamically trivial, we restrict our treatment to dynamical states with vanishing total momentum. 
We consider linear gauge conditions only. 

The Lagrangian is invariant under time-independent transformations of the Euclidean group, 

r' a = U{6)r a + u, (34) 

with U an orthogonal matrix. We define the covariant derivative D t r a = r a — £z A r a — p. Under time- 
dependent transformations r a transforms as in (|34l) and, 

e = Z-e, p' = U(e)p + u-(£-9)zAu, (D t r a )' = U(9)D t r a . (35) 

Substituting r a by D t r a in we obtain a Lagrangian that is quasi-invariant under the transformations 
(PU, and which has the form C + £ p , with C given by © and, 



^ N N N 

C p = - 2J m a p 2 - p- } j m a r a +£2-}^ m a r a A p. (36) 



a=l 



If we choose the gauge conditions £ = = p we recover the Lagrangian JQ), constrained by the eqs. of motion 
for £ and p in this gauge, 



N N 



^2m a r a Ar a = £ z , m a r a = 0. (37) 



a=l 



These constraints are first class. In the quantum theory they restrict the state space of the theory, L z ip = £ z 4>, 
^2„=i = 0, analogously to Gauss law (|A. 8|) in QED. 

We can now proceed along the same lines as in section |3 imposing on the system the gauge conditions, 



&({R a })=0, C({i? a })Ei^ m ^ = 0, (38) 

(3=1 

with & defined in l|6a|) and M = X^=i m a- ISHt defines a reference frame in a particular state of rotation, 
with origin at the center of mass. Like in section[3]we denote vectors referred to this frame by capital letters. 
The gauge conditions l|38|l are not mutually consistent unless 6 is translation invariant, 

N N 

m a A a = = m a B a . (39) 

a— 1 a— 1 

From the equations of motion for p and £ in this gauge we obtain p = and the expression 10 for £. The 
momenta conjugate to R a are given by (|10fl in terms of velocities. They satisfy and also (£({i7 Q /m Q }) = 
1 JT 7 = 0. The classical Hamiltonian is given by Q12[l. 
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The non-vanishing quantum commutators among coordinates and momenta in this gauge are, 



[x p , n Xy ] =i[6 01 - - , [y,,, n Yi ] = i [b M - -Z^l ~ ^ 



(40) 



A realization of this algebra in terms of first order differential operators with constant coefficients can be 
obtained as in (|l(j[l . 

Id lm a ^ d 1 m a A a ^ / 9 9 

P P (41) 

N N 

Id 1 m„ \ -* d 1 m n B n \ - ( . d „ d 



These operators satisfy the constraints <3({ JT a /m Q }) = = C({JT a /m a }). The realization of the residual 
angular momentum as a differential operator is the same one as in section [2 eq. 117[) . since the extra terms 
in (|41J) with respect to (|16|l do not contribute to A = ^2 a= i(X a IlY a ~ Y a IIx a ) due to the gauge condition 
£ = 0. The Hamiltonian operator has the same form as in (jTHJ), but now with the momentum operators Tl a 
from H41fl . The inner product obtained with the Faddeev-Popov procedure is, 

N 

w#>= / n rf2i? /3 <5 ( e ) ,s(2) ( c ) e ( Q w({- R «})^({^})- ( 42 ) 

3 0=1 

The quantum mechanical potential in this case is the same as in 



6 Quasirigid systems and the Eckart frame 

We assume now that the potential energy V (with U = 0) has a minimum for some configuration {z a } of 
the system, such that z a ^ z@ for some a ^ j3, and that Vo = V({2 Q }) < V({r 7 }) for all configurations 
{r 7 }. Due to the invariance of V under the Euclidean group E% any configuration {z' a } related to {z a } 
by a transformation of the form Ij34(l is also a minimum. Denoting by Ai\> the manifold of configuration 
space defined by V({r 7 }) = Vo, we assume that the quotient Mv/E 2 is a discrete set. The configurations 
of minimal potential energy are therefore rigid. In this section we discuss the quantization of the small 
oscillations of the system about these rigid equilibrium configurations. We will denote by {Z a } the unique 
(up to discrete degeneracy) minimum of V satisfying, 

N N 

^ m a Z ax Z ay = 0, ^ m a Z a = 0. (43) 

a— 1 a— 1 

The small oscillations of the system are therefore described by trajectories of the form, 

r a (t) = z a (t) + 6r a (t) with z a (t) = U{t)Z a + u (44) 

for some orthogonal matrix U(t) and u appropriately chosen so that 8r a (t) are small with respect to 
their characteristic scale for all t. We do not assume, however, that the velocities Sr a are small. Since 
we restrict ourselves to states with vanishing total momentum, the translation vector u in (|44|) must be 
time-independent. 

It is convenient to apply the inverse of the gauge transformation defined by the second equation in (|44Jl 
in order to switch to a reference frame, the "body frame" of the rigid equilibrium configuration, so that 

r a (t) = Z a + Sr a (t). (45) 
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This fixes the gauge only to leading order in 8r a . We fix the residual gauge freedom by imposing a gauge 
condition on Sr ai which amounts to correcting the definition (|45|l of the reference frame by small quantities 
of first order. We choose the origin of the reference frame at the center of mass, so to first order in 6r a the 
gauge conditions must be of the form (|38fl . The choice of the coefficients A a , B a is arbitrary as long as (|39fl 
is satisfied. We then have, 

R a (t) = Z a + SR a (t), &({SR a }) = 0, €({8R a }) = 0. (46) 

The instantaneous principal axes frame of sectional for example, is defined to first order in SR a by setting 
j4q: = and B a == Z ax in 146f) . Translational invariance of 6 is ensured by the last equation in l|43|) . 
Furthermore, not all A a and B a can vanish in this case, since we assumed Z a ^ Zp for some a ^ [3. 

The residual angular momentum in a general linear gauge (|46|l is given in the classical theory to first 
order in SR a , by 



N N / N N 

A = z ■ ^2 m aZ a A SR a + z ■ ^2 m a 5R a A 5R a - £ I V\ 2 + 2 ^ m a Z a ■ SR a 



a=l 



m 2 

with y\ 2 = Y^p=i m ()Z 2 () an d 



t m l a}) [mZ a }) + 2Q({SR a }) 



(47) 



£ = =]v— -J ( 2 ' £ m <* R <* A ^«- £ *) = ip( 1 -iiS ' ^ 

/AT \ 1 W 

X I £ ^ m Q Z Q A <5fi Q -41+ • 51 m «^a A SR a + 0(SR 2 a ). 

\ a— 1 / a — 1 



(48) 



The expression for A in the linearized principal axes gauge is obtained by substituting A a = Z ay and 
B a = Z ax in l|47l) . In general, A does not vanish at the equilibrium positions SR a = 0. In order to make A 
of first order in SR a we impose instead the gauge condition 



N N 



m a Z a A R a = ^2 m aZ a A SR a = 0, (49) 



a=l 



which corresponds to (|46(l with A a — —Z ay and B a = Z ax . The gauge condition Ij49(l defines the Eckart 
frame ^2j- The general expression l|47|) for A, in this gauge simplifies to, 

N 

A = z ■'^2m a 5R a A8R a . (50) 

a=l 

Fixing a gauge in which A is of first order in SR a as in (|50|l is a necessary condition to satisfying Casimir's 
criterion for the decoupling of rotational and vibrational degrees of freedom in a quasirigid system in low 
orders in perturbation theory |12|. As shown above, such condition is not fulfilled by the principal axes 
frame \1'2\ . 

In Eckart gauge the momentum operators II a conjugate to SR a are given by <|41[) . with the values of 
A a , B a corresponding to H49J1 . and with derivatives d/dX a , d/dY a substituted by d/dSX a and d/dSY a . 
The operator A is then of the form, 

A = - L [ sx °mz - 5Ya dsx-J ~ i w \ { Zxa ds?: ZYa dsx-J ■ (51) 

a— 1 v 7 a— 1 v / 

Its coefficients are of first order in SR a . The Hamiltonian operator is obtained by substituting these 
expressions for n a and A, into JTHJ, with Q({R a }) = 0.{{Z a }) + Q({SR a }) = V\ 2 + £^ =1 m i Z i ' SRy. 
The inner product, finally, is given by 142|) . 
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6.1 The case N = 2 



As a very minimal verification of the formalism we consider a system of two particles interacting through 
an elastic potential that models a spring with rest length a, V = fc/2(|ri — r 2 \ — a) 2 . We verify that the 
Hamiltonian operator of the previous section, with the corresponding inner product, leads to the correct 
energy spectrum and wave functions in a semiclassical expansion for large a. Since the residual angular 
momentum A vanishes in a translationally invariant system for N < 3, this simple example does not provide 
an illustration of the role of terms linear in l z in 7i, nor of Casimir's condition. In this section we restore h 
in all expressions. 

We choose the minimum Z\^ 2 of the potential as, {Z\ i2 ) x = iarr^i/M, {Z\^)y — 0, with M = m\ + m 2 . 
These Z\^ 2 satisfy 1)43(1 . The gauge conditions defining the Eckart frame are then, 

&({SR a }) = a f i(SY 1 - SY 2 ) = 0, £({5R a }) = ^SRi + ^SR 2 = 0, (52) 

M M 

where ji is the reduced mass. Together, (|52|l imply SY\ = = SY 2 . The Faddeev-Popov determinant £2 and 
9\ 2 are given by, 

= fia(a + SX 1 - 5X 2 ), m 2 =fia 2 . (53) 
The momentum operators (I4H are then, 

„ h ( m 2 d mi d \ „ „ „ , 

n ^ = l{Wdsxril9Sx-J = - IIx - n «= = n *> ^ 

consistent with the gauge conditions ((52(1 . In terms of these operators we write the Hamiltonian 1(22(1 as, 

* - 5sr^. + 5=** + ^<»- - *> 2 + w (« - i) ■ < 55 > 

where k = [no 2 , and the last term gathers the centrifugal and quantum potentials. Using IIx 1 = — I7x 2 , we 
rewrite 7i as, 



'H - Wn + H t . h( u = ±-II 2 Xi + ^-(*-Y, - YYj r 



fr 2 SK 2 / 2 1\ Tkj 2 ^2 1\ ( A n IJw, 



(56) 



In 7ii we denoted e = y h/ (fiuja 2 ), which is our perturbation expansion parameter. Taking into account 
((52(1 . 1(53(1 and the fact that the Jacobian has been absorbed in the wave functions, the inner product (14211 
takes the form, 



(M) = — // d5X l dSX 2 S[^5X 1 + -?-6X 2 )Q(a + 6X l -5X 2 )(ftiP)(X 1 ,0,X 2 ,0). (57) 



1 /*/* r /fill „, 7712 

Using 1(56(1 and 1(57(1 we can compute the perturbative expansion. Since \m\5Xx + m 2 5X 2 ,IIx 1 ] = 0, the 
eigenfunctions of Ho depend on 5X 12 only through 5X 1 — SX 2 . The eigenvalues Em) n and eigenfuntions 
<fr( ) n (5Xi — SX 2 ) of Ti.Q are then those of a one-dimensional harmonic oscillator. 
To order C(e 2 ), the perturbed energies can be read off the expression l(57?|l for Tii, 

fit*) 2 ' «•> 1 



E„ - E {0)n + E(i) n , E( X )„ = — e \£ e - - I . (58) 
Using the inner product 1(5 7|). the perturbed wave functions are found to be, to C(e 3 ), 



(/> n (SX 1 - 8X 2 ) = <£(o)„ + y ^ - yj^^^n-i - y^y^0(o)n+ij > (59) 

where on the r.h.s. we omited the argument (SXi — 5X 2 ) of wave functions for brevity. These results, ((58(1 
and 1(59(1 . agree with a conventional perturbative calculation as they should. 
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7 Final remarks 



We considered above the classical and quantum dynamics of two-dimensional many-body systems in rotating 
reference frames, in a gauge invariant approach ^ [S] ■ Our treatment parallels the formulation of gauge field 
theories in non-covariant gauges in the Schrodinger representation EI E| > an d generalizes to TV-body 
systems the analysis of the one-particle case of 0] ED] ■ 

The gauge-invariant approach allows us to deal with constrained degrees of freedom without necessarily 
solving the constraints. Yet, that approach entails also a reduction of configuration space by the elimination 
of the angular degree of freedom. The remaining, constrained dynamical variables only span the reduced 
configuration space. Our formalism has a direct physical and geometrical interpretation which we have 
tried to emphasize. In the gauges of sections I3HH1 it 1S n °t difficult to obtain the metric tensor on the gauge 
hypersurface in terms of coordinates ({i? Q }, 9) (and p if there is translation invariance). In those coordinates 
the metric has as many zero eigenvalues as the dimension of the gauge group, its kernel comprising the 
subspace orthogonal to the gauge surface. From the restriction of that metric tensor to the gauge surface 
the kinetic energy operator in the quantum theory can be constructed by applying the usual expression for 
the Laplacian in curvilinear coordinates. 

In sections we gave a derivation of the quantum theory in rotating frames. A detailed discussion is 
provided there of Gribov ambiguities and of the commutator algebra, which are essential for the obtention of 
the inner product and of the momentum and Hamiltonian operators. Those issues seem to us to have been 
neglected in the previous literature. Also discussed in detail is the residual angular momentum operator 
A. In linear gauges its eigenvalues are integer. In the quadratic gauge of section 0] its eigenvalues depend 
on the dynamical variables {R a } though, remarkably, only through Q and R 2 (see (|B.13(1 . (IB . 81 and Q23J1). 
which are measures of the instantaneous shape of the system. The relation between wave functions in the 
laboratory frame and in rotating ones, which we treat rather briefly, is summarized in the expression for the 
inner product in the different gauges, JSJ), |2"T]l. i|3*3*)l and lf4"2T) . Such relation is best understood by going 
through the succesive steps of the derivation of the inner product by the Faddeev-Popov technique. For 
systems of identical particles the formalism is clearly symmetric under permutations of particles if the gauge 
conditions are chosen so that they are symmetric in the position vectors {R a }. 

In the case of quasi-rigid systems we recovered in section El some of Eckart's classic results [Ej- In 
particular, in the classical theory the residual angular momentum A becomes manifestly small of 0(6R) in 
Eckart gauge. Thus, in a perturbative expansion the term L Z A in Ti. is of higher order than the term L 2 . 
This leads to a decoupling of rotational (L z ) and vibrational ({SR a }) degrees of freedom in low orders in 
perturbation theory. (In the elementary example of section T6. II in which A = 0, the decoupling is apparent 
through 0(e 2 ) in l|56[l.) Such decoupling would not be manifest in other gauges in which A ~ 0(1) |12l I21j . 
since in that case the terms L 2 and L Z A would be of the same perturbative order. That is the case of the 
instantaneous principal axes frame |12l l2~T] . 

We do not agree with the point of view of |5] (see also 0101) that, by quantizing a many-body system in 
the instantaneous principal axes frame (see eq. (II. 4) in |§]), it is possible to separate the collective rotation 
from the intrinsic dynamics. If that statement were generally valid, in the quasi-rigid case it should imply the 
decoupling of the total angular momentum from the vibrational degrees of freedom in perturbation theory. 
But, as shown in |12| . that decoupling is not manifest in the principal axes frame. 
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A Electrodynamics 

The methods used in the preceding sections are common to all (Abelian) gauge theories and, therefore, have 
a direct counterpart in electrodynamics. In this appendix we briefly discuss those aspects of electrodynamics, 
following the treatment of In terms of the field-strength tensor F^ v — d^A v ~d v A^, the Lagrangian 

density of the e.m. field coupled to an external current density has the familiar expression, 

C = --F^F^-fA ll . (A.l) 

We assume that the current is conserved, d^jp = 0, so that the action of the system is invariant under 
the e.m. gauge transformations A^' = A^ — d^A, with A = A(x) an arbitrary function of the space-time 
coordinates. We consider only fields, currents, and gauge-parameter functions vanishing at spatial infinity. 
Furthermore, we assume that is such that there is an inertial frame in which d j — 0, and therefore 
V • j = 0. We choose that reference frame to formulate the Hamiltonian formalism. The extension to the 
general case d°j° ^ 0, which involves time-dependent constraints, is not necessary for our purposes. 
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Weyl gauge Since the time derivative of A does not enter £, it is an auxiliary field. By means of a gauge 
transformation any field configuration can be brought into the form A = 0. We denote fields in Weyl gauge 
by V fl (x), with V° = 0. In this gauge the Lagrangian density reduces to, 

C=\v 2 -\{VAVf+j-V. (A.2) 

C must be supplemented by the equation of motion for V° derived from l|A.l|) . which is Gauss law, 

V • V = -f, (A.3) 

with V = —E. In the Hamiltonian formulation eq. I|A.3|) is written in terms of the momenta canonically 
conjugate to V, ~E, and constitutes a primary first-class constraint. It does not give rise to further 
secondary constraints. The Hamiltonian density is given by, 

U = ^E 2 + ^{V AV) 2 -j-V, V-E = f. (A.4) 
The Poisson brackets are canonical, 

[&(t,x),V k (t,y)] p = P k S(x - y). (A.5) 

All other brackets among basic dynamical variables vanish. The theory in this gauge is invariant under 
gauge transformations with a time-independent parameter A(x). The canonical generator of that symmetry 
is the l.h.s. of Gauss law, V • E — j (or just V • E). We have the Poisson brackets, 

[V • E(t, x) - f{x), H] p =Q, [V • E(t, x) - f{x), V k {t, y)) p = d k S(x - y), (A.6) 

where H = J d 3 xH(t, x). The first eq. in l|A.6|l expresses the consistency of the constraint with the dynamics, 
and the second shows that the constraint is the infinitesimal generator of the residual gauge symmetry. 

In the quantum theory in the Schrodinger representation the field operators are time- independent, the 
basic commutators are obtained from l|A.5(l . and the Hamiltonian density operator is given by, 

«(«) - -I wh^VHx) + A v{x))2 - ^ ' v{x) - (AJ) 

Gauss law is imposed as a constraint on the state space of the system, 

i v -w^)^ v] - f{xmn (A - 8) 

where ^[V] is the wave functional. 

Coulomb gauge Any field in Weyl gauge can be transformed into Coulomb gauge, V • A = 0, and 
conversely, 

Af i (t,x) = V fi (t,x) + d> 1 A(t,x), A(t,x) = ^ T (V-V(t,x)) = -— / d 3 x' V ■ X \ (A.9a) 

V 47T J I x' - x\ 

V' M (t,x) = A> l (t,x)-di M \(t,x), \(t,x)= [ dt'A {t' lX ). (A.9b) 



j t 

A in (|A.1|I is an auxiliary field which is determined by its Lagrangian equation, 

1 
V 

Once the gauge has been fixed, we can substitute this expression back into the Lagrangian, to obtain, 



A = -^f. (A.10) 



C= l -A 2 -\{V AAf+ l -f^f+j-A, V-A = 0. (All) 



15 



The electric field is 



E — Et -+- Ei 



E L = -VA° = V^V E T = -A. 



V 2 



(A.12) 



The momentum conjugate to A is 77 = —Et- Substituting l|A.9b(l and (|A.12J| into l|A.4|) . we get the classical 
Hamiltonian density in this gauge, 



n 



n z 



-(V A A) 2 



2 J V 



:j° - j ■ A, V • A = = V • 77. 



(A.13) 



Furthermore, from the inverse transformation l|A.9a|) and (|A.12|I . and the Poisson brackets (|A.6|I . we can 
compute the brackets among coordinates and momenta in this gauge. The corresponding equal-time quantum 
commutators are, 

1 



A^t,x),E k T (t,x') 



-6 



V J V 'fe S(x-x' 



(A.14) 



The gauge condition V • A = and the derived relation V • 77 = constitute second-class constraints, valid 
as operator equations. Their associated Dirac brackets are given by (|A~14)l . 

A realization of the commutator algebra l|A.14|) in the Schrodinger representation can be found as in 
sections [3] and 01 



n j (x) = - 



iSA>(x) i 



V 2 8A k {xY 



(A.15) 



where the repeated index k is summed from 1 to 3. We can express the momentum operators in Weyl gauge 
in terms of the operators (|A. 1 f>|> by means of the chain rule, 



5V m (x) 



= / d A y 



,6A k (y) 



8V m (x) 5A k (y) ' 



From (|A.9a|l we see that, 

6A k {x') 
5V m (x) 

Replacing (|A.17(I into l|A.16|) . we get, 



5 km S{x - x') 



-V' fe V 



4tt 



in m {x), 



\x' — XI 



(A.16) 



(A.17) 



(A.18) 



SV m (x) 

and therefore, from the Hamiltonian density operator i|A.7|) in Weyl gauge, we obtain the operator in Coulomb 
gauge which has the same form as the classical density (IA.13|I with 77 from l|A.15|) . The Faddeev-Popov 
Jacobian in this case is field- independent, as can be seen from (|A.17ll . which is why it does not appear in 
the Hamiltonian. 



B Eigenvalues and eigenfunctions of A 

In this appendix we discuss the eigenvalues and eigenfunctions of the residual angular momentum A = 
{XpHyp — YpIIxp)- Since the definition of 77 Q depends on the gauge conditions, the form of A as a 
differential operator is different in the cases of linear and quadratic gauge conditions. 

B.l Linear gauges 

In the linear gauges of sectionsEHElandElA is given by (jl7|l . In order to find its eigenvalues and eigenfunctions 
we apply the method of characteristic lines [22] (which can also be applied, of course, to the standard angular 
momentum operator in Cartesian coordinates). It is therefore necessary to obtain first the classical orbits 
generated by A. 
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Classical orbits The orbits generated by the classical A on the gauge surface in configuration space are 
described by the equations, 

^X = [A,X 7 ]p, ^Z = [A,y>, 6({fi 7 }) = 0. (B.l) 

With the Poisson brackets given in 114|l (as quantum commutators), we obtain the solution to (|B.1|) as, 

X 7 (a) = cxp ([A, .] P ) X 1 = + cos a (x 7 - + sin a (V 7 + ^ 

F 7 (a) = exp ([A, •} P ) Y 1 = + cosa (V 7 + ~ sina ^A 7 - 

If 6({i? 7 (0)}) = 0, then 6({i2 7 (a)}) = for all a. In that case is constant along the orbits. The center 
of mass also vanishes for all a if it vanishes for a = 0, as long as S is translationally invariant, as explained 
in sectional 



Kernel We consider first functions C satisfying AC = 0. From (|B.2|) . we see that the quantities 

are constant on the classical orbits. Thus, any / = f{{pp}) satisfies A/ = 0. Furthermore, from l|15fl we 
have AQ = i&, so that for any function / = /(£}) we get A/(Q) = f'(£})& which vanishes on the gauge 
surface. We therefore have AC = if C = C({pp}, 0) and 6 = 0. 

Eigenvalues and eigenfunctions Next, we consider the eigenvalue equation iA*f? — iAty, with A ^ 0. 
The characteristic lines of this equation are defined by the differential system 

dXl </Xx dYl dYN =^r, 6({iM) = 0. (B.4) 



Writing ^ — exp(iAa), the solutions of (|B.4() are given by the classical orbits {X 1 (— a), Y 1 (— a)} from l|B.2|) . 
A solution a to the eigenvalue equation is obtained by inverting a relation of the form, 

G(Xi(Q), . . . , X N (0), Y 1 (0), . . . , y n (0)) = 0. 

A possible choice is G = 5^(0) + -sfe-- Setting Y 1 (a) = Y y and X 1 (a) = X 1 , we are led to 

m 2 Y 7 + 0A 7 
tana ^^A 7 -£HV 

From here we obtain the family of cigenfuntions, 

(n \ N 

i -^7 a 7 ) > witn X! = ^ ^ B ' 5 ) 
7=1 / 7=1 



The proportionality "constant" C in l|B.5|) can be any function belonging to the kernel of A, as described 
above. For ^ to be single- valued A 7 , and therefore A, must be integers. 

B.2 Instantaneous principal axes gauge 

In the quadratic gauge of section 01 A has the form (|3f)|> . We proceed as in the previous case, starting from 
the classical orbits generated by A. The notation is the same as in section 
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Classical orbits In this gauge we have the following equations for the classical orbits on the gauge surface, 

^ = [a,x 7 ] p =(i + 2 4)y„ ^ = [ a,y 7 ] p = -(i-^)x 1 , S({R y }) = 0, (B.6) 



Their solution is, 



?2 , o^>\ 1/2 



' R + 2Q 

X 1 (a) — exp ([A, X 1 = cos (fia) X 1 + ( — ^ — ) sin (Qa) Y y , 



(B.7) 



R 2 -2Q 

( R 2 — 2Q\ 

Y 7 {a) = exp([A,»] P ) F 7 = - ( fi2 + 2 Q j sin (^«)^7 + cos {fla)Y 7 , 
with, 

" = ('-1?f- < B -«> 

The condition S = is preserved along these orbits and, if it is satisfied, Q and R 2 are constant on them. 
Notice that we require Q > due to the Gribov ambiguity (see section 0} and, therefore, by definition 
< 2Q < R 2 . 



Kernel The quantities 



p^v-w x t + ( 1+ r¥ Y * (R9) 



are constant along the orbits l|B.?jl . Thus, if S = 0, any function C = C({p 2 }, Q, R 2 ) satisfies AC = 0. 

Eigenvalues and eigenfunctions The characteristic lines of the equation zAvP = iX^, A ^ 0, in this 
gauge are the solutions to the system, 

dXi dX N _ dY 1 dY N _ 

together with the gauge condition 5({i? 7 }) = 0. As in the previous section, we set = exp(iAa). With 
this, the solution to (|B.10|I is given by {AT 7 (— a), Y^(— a)}, see (|B.7|) . Inverting the relation Ky(0) = 0, we 
get, 

1 f{R 2 + 2Q\ 1/2 YA 

a, = -arctan [[jf^ ) ^J- 

We then have, 

i^A 7 a 7 j, A = ^A 7 . (B.12) 

7=1 / 7=1 

For 4" to be singled-valued we need, 

N 

A 7 = n 7 f2, A = ?ifi, rt = ^^n 7 , (B.13) 

7=1 

with n 7 integers. We see that the eigenvalues A are not integers in general, and that they depend through 
£1 on Q({R a }) and R 2 ({R a }), which in turn are functions of the dynamical variables {R a }. 
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